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ABSTRACT 

Aims. Despite photometry and spectroscopy of its oscillations obtained over the past 25 years, the pulsation frequency spectrum of 
the rapidly oscillating Ap (roAp) star y Equ has remained poorly understood. Better time-series photometry, combined with recent 
advances to incorporate interior magnetic field geometry into pulsational models, enable us to perform improved asteroseismology of 
this roAp star. 

Methods. We obtained 19 days of continuous high-precision photometry of y Equ with the Most (Microvariability & Oscillations 
of STars) satellite. The data were reduced with two different reduction techniques and significant frequencies were identified. Those 
frequencies were fitted by interpolating a grid of pulsation models that include dipole magnetic fields of various polar strengths. 
Results. We identify 7 frequencies in y Equ that we associate with 5 high-overtone p-modes and 1st and 2nd harmonics of the 
dominant p-mode. One of the modes and both harmonics are new discoveries for this star. Our best model solution (1.8 Mq, log 
Teff ~ 3.882; polar field strength ~ 8.1 kG) leads to unique mode identifications for these frequencies (£ = 0, 1, 2 and 4). This is 
the first purely asteroseismic fit to a grid of magnetic models. We measure amplitude and phase modulation of the primary frequency 
due to beating with a closely spaced frequency that had never been resolved. This casts doubts on theories that such modulation - 
unrelated to the rotation of the star - is due to a stochastic excitation mechanism. 

Key words, stars: chemically peculiar - stars: oscillations - stars: individual: y Equ - stars: magnetic fields - methods: data analysis 



1 . Introduction - the history of y Equ 

Rapidly oscillating Ap (roAp) stars form a class of variables 
consisting of cool magnetic Ap stars with spectral types rang- 
ing from A-F, luminosity class V, and were discovered by Kurtz 
(1982). Photometric and spectroscopic observations during the 
last 2 decades characterize roAp stars with low amplitude pul- 
sations (< 13mmag) and periods between 5 to 21 minutes. It 
is now widely accepted that the oscillations are due to high- 
overtone (« > 15) non-radial p-mode pulsations. The observed 
frequencies can be described in a first-order approximation by 
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the oblique pulsator model as described bv lKurtzl dl982l) . In this 
model the pulsation axis is aligned with the magnetic axis, which 
itself is oblique to the rotation axis of the star. It has since been 
refined a nd improved to include additi onal effects like the corio- 
lis force dBigoi & Dziembowskill2002l) . 

The oscillations in roAp stars are most probably excited 
by the k mechanism acting i n the hydrogen ionization zone 
dDziembowski & Goodd fl99q) . while the prime candidate for 
the selection mechani sm of the observed pulsation modes is 
the magnetic field (e.g . lBigot et alJfeOOOtlCunha & Goughl2000l: 
iBalmforth efailJgOQU (Bigot & Dziembowskill2002l 120031 ISaiol 
12.005b ICun ha 2006). Similar to non-oscillating Ap (noAp) stars, 
the group of roAp stars also shows a distinctive chemical pat- 
tern in their atmospheres, in particular an overabundance of 
rare earth elements, which points to stratification and vertical 
abundance gradients. Tackling the question of what drives the 
roAp pulsation, it is important to note that apart from the sim- 
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ilar chemical composition, noAp stars show comparable rota- 
tion periods and magnetic fiel ds but higher effective tempera- 
tures dRvabchikova et aLll2004l) . Thus, their hydrogen ionization 
zone lies further out in the stellar envelope, which prevents ef- 
ficient driving of pulsation. All together, roAp stars have given 
new impulse to ast eroseismology and 2D mapping of pulsation 
dKochukhovl |2004|) . as well as detailed 3D mapping of (mag- 
netic) stellar atmospheres, has become possible. For reviews 
on primarily observation al as pects of roAp star s we refer to 
iKurtz & Martinez! d2000t) and IKurtz et al ] (12004 and on theo- 
retical aspects to lShibahashild2003l) and lCunhal(l2005l) . 

y Equ, (HD 20 1601, HR8097, A9p, m v = 4.7) was the 
6th discovery as a roA p star. It was lon g known to have a sig- 
nificant magnetic field dBabcocklll958l) which is variable with 
a per i od of ~ 72 years dBonsack & Pilachowskil Il974t IScholj 
119791: iLerov et al.lll994l) . Rotation, a solar-like magnetic cycle 
or precession of the star's rotation axis due to its binary com- 
panion are candidate mechanisms fo r this cyclic variatio n. The 
binary hypothesis was supported by IScholz et al.l dl997l) . who 
observed a temporary drop in radial velocity for y Equ dur- 
in g 4 consecutive nights. Their re sults are however contradicted 
by lMkrtichian et al.l (1 19981 Il999l) . Magnetic field data spanning 
more than 58 years indicate that the variation is most like ly due 
to rotation and therefore P mag = P m t dBvchkov et al.ll20 06T). This 
explanation is also supported by many spect roscopic observa- 
tions showing very sharp absorption lines (e.g. lKanaan & Hatzesl 
119981) which are typical for a slow rotator (or a star seen pole- 
on)^ 

iKurtzl (fl983h was the first to detect a pulsation period 
of 12.5 min (1.339mHz) with an amplitude varying between 
0.32 and 1.43 mmag and speculated about a rotation period 
of 38 days, inconsistent with the magnetic field measurements 
mentioned above. Aside from rotation, beating with a closely 
spaced frequency has been proposed as a cause, for which the 
present paper gives further evidence. The p u lsation was con- 
firmed shortly after Kurtz's study by IWeissI dl983l) . iBvchkovl 
dl987l) reported first evidence for radial velocity ( RV) varia- 
tions, but it took two more years for a clear detection dLibbrechj 
1 1988b . Although unable to detect the 1.339 mHz oscillations 
reported by Kurtz (1983), Libbrecht discovered three frequen- 
cies at 1.365 mHz, 1.369 mHz and 1.427 mHz. He suggested 
that the amplitude modulation observed in the spectra of roAp 
stars may not be due to closely spaced frequencies, but rather 
caused by short mode lifetimes in the order of ~ 1 d. He con- 
cluded that both peaks therefore b elong to a single p-mode 
oscillation. The follow-up study by IWeiss & Schneider! dl989l) 
aimed at performing simultaneous spectroscopic and photomet- 
ric studies but failed to confirm y Equ's pulsation. The latter 
authors applied a correlation technique using more than 100 A 
wide spectral ranges in order to improve the S/N ratio relative to 
previous techniques, based on individual lines. This approach 
turned out to be inapplicable, considering the peculiar abun- 
dance s tructure of roAp star atm ospheres which is now well 
known. iMatthews & Scottl (1 19951) then claimed to have detected 
variability with rather large radial velocity amplitudes but their 
results were inconsistent with previous findings. 

Detailed photom etric observations were conducted by 
iMartinez et al.l d 19961) using a multi-site campaign in 1992, span- 
ning a total of 26 nights. Their results also suggested limited 
life times of pulsation modes, because additional frequencies ap- 
peared in their analysis of individual nights. However, all three 
frequencies detected so far were confirmed, and prewhitening 
of 1.366 mHz revealed a fourth eigenfrequency at 1.397 mHz. 
Using these four frequencies they concluded a spacing of con- 



secutive radial overtones of Av ~ 30 ^/Hz based on the asymp- 
totic theory fo r low-degree, high-overtone p-mode pulsations 
(lTassoull ll990). They did not detect any evidence for non-linear 
behavior of the eigenfrequencies (i.e. harmonics) which, as will 
be shown in this paper, is indeed present in the Most data. So far 
their observing campaign was the most recent photometric in- 
vestigation of y Equ. A period of spectroscopic studies focusing 
on pulsation, abundan c e and stratification analyses followed. 

iKanaan & Hatzesl (1 19981) reported on RV amplitudes for 
chromium and titanium and speculated that this is due to 
their concentration close to the magnetic - hence also pul- 
sation - poles, while other elements (e.g. iron) showing no 
RV variations may be con centrated at the magnetic equator. 
iMalanushenko et all (1 19981) independently arrived at a similar 
conclusion and they are the first who identified the rare earth el- 
ements (REE) Pr in and Nd m to show the largest variations (up 
to 800 ms -1 ). Most other atomic species had very low or non- 
measurable RV variations and the authors concluded that previ- 
ous spectroscopic observations failed to detect significant radial 
velocity variations due to the chosen spe ctral regions and/or low 
spectral resolution. ISavanov et al.l d 19991) further clarified the sit- 
uation by drawing attention to a possible inc orrect identification 
of spectral lines in IKanaan & Hatzesl (fl998 ) and concluded that 
Pr m and Nd m were responsible for the large RV amplitudes dis- 
cussed in their analysis. 

Further investigations bv | Kochukhov & Rvabchikoval(l200ll) 
and lRvabchikova et al.1 (|2002) impressively explained amplitude 
modulations of different elements and even ions of the same el- 
ement by lines being formed at different atmospheric depths. 
The authors also tried a first mode identification based on rather 
short data sets of line profile variations (LPV s) and they argue 
for I - 2 or 3, and m - —6 or —t + 1 modes. Shib ahashi et al.l 
(|2004|) disagreed with this analysis and suggested a shock wave 
causing the observed L PVs which in turn was questioned by 
iKochukhov et al.l d2007l) who argued that the pulsational velocity 
should not exceed the local sound speed. The latter authors pro- 
pose a modified oblique pulsator model where LPVs are caused 
by pulsation velocity fields superposed by sinusoidal line width 
changes due to convection, but which previously was thought to 
be suppressed. In any case, they agree with Shibahashi et al. that 
the identification of the primary frequency as an I = 1, m = 
mode is still the most likely explanation, which is also sup- 
ported by this paper. T he issue is still being debated though 
ilShibahashi et alj|2007l) . 

A still unsolved is sue are magnetic field v ariations synchro- 
nized with pulsation. iLeone & Kurtzl (120031) showed evidence 
for this, but subsequ e nt investigations (IKochukhov etail 120041: 
iBvchkov et al.l 20051: Savanov et al.l 12006 ) could not confirm 
their findings. iHubrig et al.1 (120041) corroborate this null detec- 
tion for y Equ with ESO-FORS1 data. 

Despite a history of more than 25 years of observations, the 
pulsation frequency spectrum of y Equ remains poorly under- 
stood. Several frequencies have been published, but up to now 
these have never been observed simultaneously. This is where 
Most steps in, providing the most complete high precision pho- 
tometric campaign ever conducted for y Equ and covering con- 
tinuously a timespan of 19 days. 

2. Observations and data reduction 

2.1. Most observations 

Most is a Canadian space satellite designed for the detection 
of stellar variability with amplitudes down to several ppm on 
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Table 1. Reduction statistics for y Equ. A total of 49213 frames 
were gathered. ACS: Most Attitude Control System. 



Reduction step 


# rejected 


% rejected 


saturated pixels 


1 


0.002 


saturated ACS errors 


1637 


3.33 


deviating image geometry 


21679 


44.35 


excessive number of cosmics 


1002 


3.87 



time scales up to several days (I Walker et al.ll2003t) . The satel- 
lite was launched in June 2003 and has since proven to deliver 
photometric time series of unprecedented precision. It carries a 
Maksutov telescope with an aperture of 15 cm, uses a broadband 
filter (350-700 nm) and operates in a low Earth orbit following 
the terminator in an almost polar orbit. Most can observe stars 
in the continuous viewing zone for up to 6 weeks with a point- 
ing precision of ~ 1". An array of Fabry lenses is used to sup- 
press effects of satellite jitter in the data of primary targets of 
V < 6 mag. y Equ was observed for 19 days from 07/28/2004 
to 08/16/2004 with exposures (integration time =11 sec) taken 
continuously once every 30 sec. As y Equ was one of the first pri- 
mary targets observed by Most, stray light, produced by earth- 
shine, and other instrumental effects were still being investigated 
at this time with the aim of improving the observing procedures 
for future runs. In fact, this data set was the basis for the develop- 
ment of the C reduc tion pipeline for Most Fabry Imaging targets 
dReegen et al.ll2006l) . 

2.2. Data reduction using decorrelation 

Stray light effects are corrected by computing a correlation be- 
tween target and background pixels distinguished by a fixed 
aperture, where target denotes pixels inside the aperture and 
background characterizes pixels outside the aperture, which are 
assumed to contain no stellar signal at all. A "neutral" area flag, 
for pixels not used for either purpose, is meant to exclude pixels 
at the transition zone close to the edge of the aperture. 

As it turned out throughout the Most mission, each target 
star's photometry suffered specific problems which needed to be 
addressed individually by optimizing the reduction parameters. 
Due to bad data quality at the beginning of the y Equ-run about 
1060 frames had to be rejected. Additional frames needed to be 
eliminated due to pointing problems in the early days of Most 
operations which resulted in distorted Fabry image geometries 
identified by comparison to a mean normalized Fabry image. If 
the pixel values in an individual (normalized) frame deviate by 
more than g cr, where g is an integer, from the pixel values of the 
mean normalized reference image, the frame is rejected. Since 
the mean reference image changes after the elimination process, 
this step is repeated iteratively until no more frames need to be 
rejected. Statistics of the reduction process are shown in TableQ] 

Tests showed that the number of images in which pixels are 
deviating by more than 5 cr first decreased as expected, but in- 
creased after about ten iterations (see Fig.Q] upper panel) and 
later decreased again, indicating two separate causes for devi- 
ating Fabry image geometry. Fig. [2] shows the raw light curve 
after correcting for cosmic rays, phased with the orbital period, 
clearly indicating which phases are responsible for the image ge- 
ometry deviations. Not surprisingly, the majority of these rejec- 
tions concern exposures taken during high stray light phases and 
therefore cause regular gaps in the final data set (see Fig.[3j. To 
ensure that this reduced duty cycle does not impair the identifi- 
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Fig. 1. Upper panel: Iterative elimination of frames with devi- 
ating image geometry. The increasing number of rejections with 
higher iterations followed by a second decline indicates two sep- 
arate mechanisms responsible for distorting the Fabry image ge- 
ometry. 

Lower panel : Window function of the data sets including (solid 
line) and excluding (dotted line) frames with deviating image 
geometry. The dotted graph has been shifted by 0. 1 for better 
visibility. 



cation of intrinsic frequencies, the analysis was repeated with the 
data set obtained after the first iteration step in the image geom- 
etry evaluation, which eliminated only 807 frames. In the lower 
panel of Fig.[T]we compare the spectral windows of both reduc- 
tions, showing that the rejections mainly affected orbit-induced 
artifacts (/orbit = 14.19d _1 ). The final analysis did not yield a 
different set of intrinsic, but considerably fewer instrumental fre- 
quencies. 

2.3. Data reduction using "doughnut" fitting 

Another (less invasive) reduction method was developed to ver- 
ify the consistency of our frequency analysis. In a first step an 
average Fabry frame is determined from all frames obtained dur- 
ing orbital phases with very low stray light signal. The mean 
intensity, derived from pixels which are defined as background 
("sky") pixels, is subtracted from all pixels. The resulting frame 
is scaled to the mean intensity of the N pixels with the high- 
est intensity values, where N is of the order of 100. The re- 
sulting frame serves as sort of normalized point spread function 
(PSF) for a Fabry image (Most's Fabry lenses produce images 
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0.4 0.6 
orbital phase 

Fig. 2. Upper panel : The raw light curve after cosmic ray cor- 
rection, phased with the Most orbital period. The stray light con- 
tamination is mostly confined to orbital phases - 0.4. 
Lower panel : Same as in the upper panel, but after correcting 
for deviating image geometry. 




1680 
HJD - 2451545 



Fig. 3. Final light curve using decorrelation resulting in a duty 
cycle of ~ 53%. All long-term variations with periods > 1 d 
have been prewhitened. The inset shows a subset spanning about 
50min together with a 10 point (~4min) running average. The 
~ 12min oscillation is clearly visible. 



of the entrance aperture of the optics which resemble a dough- 
nut - hence the name "doughnut fitting"). As a next step a linear 
regression between pixel intensities of an image and the corre- 
sponding pixel intensities of the mean PSF frame is performed. 
For each of these images the slope k of the linear fit gives a 
scaling factor, while the offset d corresponds to the image mean 
background intensity. Each frame is replaced by the PSF frame 
multiplied with the corresponding scaling factor k and the mean 
background intensity d is added. The average intensities of all 
target and background pixels result in the target and background 
light curve, respectively. Finally, a linear fit between the back- 
ground and target light curve is determined and subtracted from 
the target light curve resulting in the final light curve. In the case 
of y Equ the latter consists of 48955 datapoints. 

The advantage of this method is the insensitivity of the linear 
regression to extreme pixel values, for instance due to cosmic ray 



hits or local stray light effects. Also, contrary to th e method de- 
scribed in section [2~2l (see also lReegen et al.ll2006l) . the data set 
is not split into subsets and processed individually, which could 
distort the low frequency signal. The entire data set is reduced 
en-block, no subsets are created which qualifies this procedure 
as a perfect cross-check for detecting artifacts due to data reduc- 
tion. The disadvantage of this method is the stray light correction 
which is less efficient as for the decorrelation method. 



3. Frequency analysis 

Time series resulting from both redu ctions outlined in section [2] 
were analyzed individually. SigSpec ( lReegenll2007l) was used to 
obtain frequencies by means of the spectral significance between 
and 360 d -1 in a prewhitening sequence down to a significance 
level of 5.46, corresponding roughly to an amplitude S/N ratio 
of 4. Only frequencies without a counterpart in the sky back- 
ground were consid ered, following the procedure described by 
iReegen et alj (120081) . Most background light curves typically 
contain several dozen frequencies between and 360 d _1 , pro- 
duced by stray light and instrumental effects. These are, in most 
cases well confined to the orbital frequency (/01-bit = 14.19d _1 ) 
and harmonics of it, to ld~' side lobes due to the passages of 
Most above nearly the same ground pattern after 14 orbits, or 
are caused by temperature drifts of electronic boards. In the 
case of y Equ we used a method similar to what is described in 
iGruberbauer et al.l d2007l) to identify non-intrinsic frequencies. 
Finally, the remaining frequencies deduced from both reduction 
methods were compared and, again, only matching frequencies 
were taken into account for our final analysis (Fig. [4]). 

For the early Most runs the low-frequency region is known 
to be affected by an instrumental effect producin g spectral fea- 
tures centered on 3.16d~' and multiples of it dReegen et al.l 
2006). This instrumental effect is much weaker for the back- 
ground pixels which explains why the corresponding frequencies 
remain even after comparisons with the background frequency 
spectrum and what justifies the rejection of all frequencies be- 
low 10d _1 . Some of the power excess in the low frequency re- 
gion might be still of stellar origin, but the current data is incon- 
clusive. We also found the frequencies between 10 and 35 d -1 to 
be multiples of 3. 16 d , but no other Most data set shows these 
instrumental effects at such high frequencies. Because we prefer 
a critical approach, in the end, only the six higher frequencies 
with / > 100 d -1 were considered to be intrinsic to y Equ be- 
yond doubt. 

Table|2]lists the result of our conservative analysis. Two ad- 
ditional significant frequencies, not included in the table, can be 
found very close to f\ with amplitud es of ~ 40 ppm, but in clud- 
ing them in a Period04 multi-sine fit dLenz & Bregerl2005l) fails, 
because the solution does not converge. Their nature is discussed 
in section l6\2l 

No obvious spacing, as predicted by the asymptotic theory of 
non-radial pulsation, can be found in our final set of frequencies. 
This suggests that the regular spacing of modes with the same 
degree t is disturbed, e.g. by the magnetic field, or that modes 
of different degree are excited. We therefore turned to fit our 
observed frequencies to the latest generation of roAp-pulsation 
models, but using only f\ to f$ for the fitting process, as fa and 
fi correspond to the harmonics of f\ . 

4. Pulsation models 

To compare theoretical frequencies with observed ones of y Equ, 
we have computed main-sequence evolutionary models for a 
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Table 2. The list of frequencies considered to be intrinsic to 
y Equ. g-f\ 1 cr-uncer tainty in frequency, derived according to 
iKallinger et al] d2007l) ; cr a : uncertainty in a mplitude; sig: spec- 
tral significance according to lReegen] d2007b ; 6: phase in radians 
corresponding to cos(2nfjt - 0). 
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Fig. 4. The amplitude spectrum of the final light curve (gray) 
presented in Fig.[3]and all frequencies found by SigSpec in both 
reduction methods described in this paper (black). We consider 
only frequencies with / > 50 d -1 to be unaffected by the instru- 
ment, y Equ's pulsation is obvious in the range of 100 to 125 d~' 
and the first of two harmonics of f\ is evident at ~ 235 d -1 . 



mass range of 1 .75 M - 1 .85 M with heavy element abundance 
of 0.0 15 < Z < 0.025 and OPAL opacities (llglesias & Rogers! 
119961) . A list of the individual models, and how we organised 
them into grids, is presented in Table[3] For standard models, 
envelope convection is suppressed (a = 0.0), and to mimic de- 
pletion of helium abundance in the outermost layers of a star, the 
helium abundance Y is given as Y — 0.01 + 0.27(jt2 + X3), where 
X2 and X3 are fractions of singly and doubly ionized helium, re- 
spectively (cf lBalmforfh et al1l2001l) . Nonetheless, the effects of 
convection were also tested by computing all grids except Grid 
4 with = 1.5 and homogeneous He abundance in the outer 
layers, but otherwise the same parameters. Grid 4 included a ho- 
mogeneous distribution (mixing) of He in the envelope but no 
convection. 

Nonadiabatic frequencies of axisymmetric (m = 0) high or- 
der p-modes under the presence of a dipole magnetic fie ld were 
calculated using the method described in ISaiol (120051) . Outer 
boundary conditions are imposed at an optical depth of 10 -3 . A 
reflective mechanical condition (dp/p — > constant) is adopted. 



All frequencies presented in this paper are less than the acoustic 
critical frequencies of corresponding models. 

Since the latitudinal dependence of amplitude under the 
presence of a magnetic field cannot be expressed by a single 
Legendre function, we expand it into a truncated series of com- 
ponents proportional to Legendre functions Pi with lj = 2j — 1 
for odd modes and lj — 2j for even modes, where we have in- 
cluded twelve terms; i.e., j — 1, 2, . . . , 12. The latitudinal degree 
I is not a definite quantity for a pulsation mode any more, be- 
cause pulsation energy is distributed among twelve components 
associated with P\ . For convenience we still use t represent- 
ing the lj value of the component associated with the largest ki- 
netic energy. Sometimes the distribution of the kinetic energy 
among the components is broad so that the identification of ( 
is ambiguous, and the value of C may change as the strength 
of the magnetic field changes. Also, for a mode associated with 
(, the latitudinal dependence of amplitude on the stellar surface 
may be considerably different from th at of (cos 9) as shown in 
ISaio & Gautschyl(|2004 : ISaiol(l2005h (see also below). 



5. Model fitting 

Pulsation models were calculated taking into account a magnetic 
dipole with polar field strength, Bp, up to 12 kG. These models 
were tried to match the currently estimated parameter space of y 
Equ, which is log r eff = 3. 882 + 0.01 1 K, logL/L Q = 1.10 + 0.03, 
and M ~ 1.74 + 0.03 M dKochukhov & Bagnuloll2006T) . Fig.0 
shows the position and parameters of all models involved. For 
the corresponding values and further details, we again refer to 
Tabled 

To increase the resolution of our model grid, the mode fre- 
quencies were interpolated linearly in Bp and (log T e s, log L/L Q ) 
for a fixed mass and chemical composition. The step width for 
the interpolation was chosen to be 0.02 kG in Bp and 0.00002 
in log r e ff. Since not all modes with the same degree i and ra- 
dial order n converge in the model calculations, the sequence of 
frequencies for a specific mode may not cover the entire chosen 
stellar fundamental parameter space (log T e g, log L/L , and M). 
Finally we used a single coordinate, Av, according to 



Av = 0.1349 • 



0.5 



(1) 



where p is the mean density of the stellar model. All model fit- 
ting was carried o ut for a given mass in the (Z ?p,Ay)-space with 
a^ 2 -test similar to lGuenfher & Brownl d2004l) . Since an a priori 
mode identification was not a reasonable option, each model's 
eigenfrequencies of spherical degrees I — to I — 4 were com- 
pared to the five observed frequencies. The 1 cr-uncertainties of 
the observed frequencies were estimated to be roughly 0.25 of 
the up per frequency error derived according to Kall inger et al.l 
d2007l) . The corresponding 1 cr-errors assigned to the individual 
model frequencies were assumed to be 0.2 pHz. 

The best fitted models, represented by a minimum x 2 , are 
listed in Table|4] Grid 3, as defined in Table|3] outperforms the 
other models by far. It is the only model grid that manages to 
produce frequencies that fit all five observed frequencies on av- 
erage to within the uncertainties (y 1 < 1), as it is illustrated by 
the upper panel of Fig. [6] As expected, including convection for 
Grid 3 destroys the fit, while for grids with large x 1 it has the op- 
posite effect. The best fit is located at the center of an extended 
patch. Its proximity to a genuine calculated model ensures that 
this good fit is not an artifact produced by the interpolation (also 
see Section|6]l. The mean parameter values for the^ 2 < 1 region 
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Table 3. List of all stellar evolution models computed for y Equ 
with polar magnetic field strengths ranging from to lOkG in 
steps of 0.2 kG, except of Grid 7 which was calculated up to a 
magnetic field strength of 12 kG. 



Grid ID 


M/M 


Z 


loa L/ Lrs 


log X e ff 


la 


1.75 


0.02 


1.0252 


3.8810 


lb 


1.75 


0.02 


1.0229 


3.8825 


2a 


1.80 


0.015 


1.1574 


3.8990 


2b 


1.80 


0.015 


1.1558 


3.9004 


2c 


1.80 


0.015 


1.1542 


3.9017 


3a 


1.80 


0.02 


1.0880 


3.8809 


3 b 


1.80 


0.02 


1.0871 


3.8818 


3c 


1.80 


0.02 


1.0861 


3.8826 


4a 


1.80 


0.02* 


1.0899 


3.8794 


4b 


1.80 


0.02* 


1.0880 


3.8810 


4c 


1.80 


0.02* 


1.0871 


3.8819 


4d 


1.80 


0.02* 


1.0861 


3.8827 


5a 


1.80 


0.025 


1.0284 


3.8701 


5 b 


1.80 


0.025 


1.0267 


3.8714 


5c 


1.80 


0.025 


1.0243 


3.8726 


6a 


1.85 


0.02 


1.1388 


3.8885 


6b 


1.85 


0.02 


1.1368 


3.8903 


6c 


1.85 


0.02 


1.1345 


3.8917 


6d 


1.85 


0.02 


1.1325 


3.8934 


7a 


1.85 


0.025 


1.0753 


3.8806 


7 b 


1.85 


0.025 


1.0734 


3.8819 


7c 


1.85 


0.025 


1.0715 


3.8832 
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Fig. 5. A schematic diagram of the calculated grids and their 
model parameters following the notation of Table[3] Since Grid 
3 and Grid 4 partially overlap, all grid points for Grid 3 are 
shown as + symbols, while Grid 4 is indicated by rectangles. 
Observationa l uncertainties for y Equ's estim ated position in the 
HR diagram dKochukhov & Bagnuld l2006l) are represented by 
the large rectangle. 

Table 4. x 1 values of the best fits in the corresponding model 
grids. Models taking convection into account are denoted by an 
index a. 



with the mentioned interpolation step size are log T e $ = 3.8818, 
logL/L = 1.0871 and B P = 8.1 kG. 

Fig. [7] shows an echelle diagram of the best fitted model to- 
gether with the observations. fa, the primary frequency, is iden- 
tified as an I — 1 mode, fa and fa are matched by consecutive 
t = 4 modes. The remaining two frequencies, fa and fa, are fit- 
ted by modes of degree t = 2 and t — 0. As a comparison, the 
lower panel of Fig. [6] shows the fitting results to Grid 6, a grid 
with considerably larger^- 2 values. Interestingly, while all grids 
except Grid 3 fail at delivering y 1 < 1 -results, the mode iden- 
tification remains stable for grids with x 1 < 4. In particular, the 
two closely spaced model frequencies in the vicinity of fa and fa 
are found to fit the same degrees I — 1 and 4 for Grids 3, 7a, 5a, 
and 6. This is not the case for 6a and la. Thus, in our grids there 
seems to be a tendency towards the mode identification of Grid 
3 for better fitted models. Also, convection is shown to radically 
influence the model frequencies. 

In Fig.|6]obvious discontinuities can be seen in the fitting re- 
sults. They are produced when the interpolation routine cannot 
find two modes of equal spherical degree t and radial order n in 
both models that are acting as sampling points. This can happen, 
as already mentioned in Section|4] when the lvalue of a mode 
formally changes, because the former associated spherical har- 
monic does no longer supply most of the kinetic energy. Also, 
the cyclic variation of the dampin g rate as a function of Bp, as 
shown in ISaio & Gautschvl (|2004), has an effect on the "avail- 
ability" of certain modes. Around the maximum of the damping 
rate, the kinetic energy is so broadly distributed among the dif- 
ferent ^-components that convergence of the model calculations 
starts to fail. This effect is clearly visible as the two large dis- 
continuities in each panel of Fig.[6] which be ars resemblance to 
what is presented in lSaio & Gautschvl d2004f) . 



Grid ID 


min(x z ) 


S P [kG] 


1 


41.0 


6.0 


la 


4.4 


8.8 


2 


6.0 


6.5 


2a 


21.7 


6.5 


3 


0.7 


8.1 


3a 


7.7 


2.8 


4 


8.0 


0.0 


5 


22.6 


0.4 


5a 


2.2 


4.6 


6 


4.6 


7.3 


6a 


4.6 


6.0 


7 


17.8 


0.0 


7a 


2.0 


3.7 



6. Discussion 

6.1. The frequencies ofy Equ 

As mentioned in the introduction, roAp stars are characterized 
by an interaction of a strong global magnetic field with velocity 
fields due to pulsation of still poorly known origin, all leading to 
a complex eigenfrequency spectrum. The latter provide the only 
directly accessible information concerning the internal structure 
of these stars. 

Data obtained by Most have allowed us to unambiguously 
identify 7 frequencies (Table|2]) above 1.3 mHz, including the 
first two harmonics of the primary frequency, which signifi- 
cantly exceed the mean noise level of ~ lOppm between 1.15 
and 1.75 mHz. The comparison with Martinez et al. yields a 
match for fa, fa, and fa. Our fa has never been detected be- 
fore, probably because fa and fa could not be resolved as in- 
dividual frequencies. This is also important for the discussion 
on amplitude modulation of the primary frequency, for which 
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Fig. 6. Upper panel: Color coded inverse x 1 values obtained 
when fitting the 5 intrinsic frequencies of y Equ to model fre- 
quencies of Grid 3. The vertically aligned symbols represent the 
calculated model positions - all other model frequencies are de- 
rived via linear interpolation. Discontinuities are due to a lack of 
converging model frequencies. Some modes of a certain degree 
and radial order are missing in the individual genuine models. 
Also, sudden changes of the frequency spacing due to the mag- 
netic field, can change the quality of the fit quite rapidly. 
Lower panel: Grid 6 shows a pattern similar to Grid 3 but deliv- 
ers larger^ 2 values. It gives, however, the same mode identifica- 
tion. Note the different scale of the abscissa! 



I— * M* 



0.044 0.046 0.048 0.05 

frequency mod 60 uHz [mHz] 

Fig. 7. Upper panel: Echelle diagram of the best fitted stellar 
model (x 1 = 0.744, Grid 3). The corresponding parameters are 
log7; ff = 3.8818, logL/L = 1.0871, and B P = 8.06kG. The 
model frequencies are represented by their estimated error bars, 
the numbers above denote the mode degree I. Black cross sym- 
bols indicate the position of the observed frequencies with error 
bars. 

Lower panel: A close-up of the region in the echelle diagram 
including fi, fi, fo, and f$ illustrating the quality of the fit more 
clearly. f$ is not matched to within the 1 cr-error, but f\ to f% 
dominate the^ 2 statistics because of their high significance. 



we refer to the next section. Our value for f$ is comparable to 
their v\ = 1.321 mHz (taken from Weiss 1983), but differs by 
~ 0.01 mHz. We assume that they have misidentified a 1-day 
alias of the real frequency but this is difficult to asses, since no 
frequency uncert ainties are kn own for this data set. Given that 
the data set from lWeissI d 19831) spans only three nights, the poor 
frequency resolution l/r o b s - 0.33 d _1 = 0.004 mHz supports 
our assumption. Thus, we can for the first time confirm and ex- 
pand the previous set of frequencies unambiguously. 

When discussing published amplitudes (including those in 
the present paper) one has to keep in mind that two differ- 
ent properties prevent a direct comparison. First, amplitude 
modulations are present which lead to different amplitudes 
for observations obtained at different beating phases. Second, 
it has been well known since the e ighties that amplitudes 
dWeiss & SchneideJ[T984l) and phases (IWeisslll986h depend on 
wavelength and the former are rapidly decreasing towards the 



red. The passband of Most is broad compared to the Stromgren 
system, frequently used for roAp star photometry, hence the ob- 
served Most amplitudes are intrinsically smaller. 

6.2. Amplitude modulation 

Amplitude modulation of y Equ's primary frequency has of- 
ten been mentioned in the literature and rotation or closely 
spaced frequencies have been proposed as possible explana- 
tions. Because of the relatively small effect and the limited ac- 
curacy of the data a specific modulation period could never be 
accurately determined. This situation, however, changed with 
Most. To distinguish between amplitude modulation of a sin- 
gle frequency and beating of two closely spaced frequencies 
one needs to discuss simultan eous phase and amplitude changes 
dBreger & Pam vatnvkh 2006). Figure[8]illustrates the correlation 
of amplitude and phase variations for f\ with a relative phase 
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Fig. 8. Upper panel: Amplitude modulation of the primary fre- 
quency, fx, as found by a least-squares fit with a fixed frequency 
to subsets of the light curve. 

Lower panel: Phase changes of f\ determined in the same way 
as the amplitude modulation. The dotted lines represent the 1 cr- 
confidence limits. The similarity of amplitude modulation and 
phase shift, as well as the minimum amplitude coinciding with 
the time of the phase shift turning point, suggest an amplitude 
modulation due to beating of two closely spaced frequencies 
with a beat frequency of ~ 0.07 d -1 (~ 0.0008 mHz). This corre- 
sponds to the beat frequency of fx and fi. 



shift of |, which is indicative for beating of a close pair of fre- 
quencies spaced by /b ea t ~ 0.07 d _I (~ 0.0008 mHz). 

We find / + /beat - fi- The best fitting pulsation model pre- 
dicts both fx and fx to be eigenfrequencies of the star, hence 
we conclude that the "closely spaced frequencies"-hypothesis is 
correct. It comes as no surprise that previous observations, espe- 
cially time-resolved spectroscopy, could not explain the modu- 
lation effect, since it is impossible to resolve fx with short time 
bases. The possibility of unresolved modes in such data should 
be taken into account when studying pulsation via residual spec- 
tra, produced by subtraction of a mean spectrum phased by only 
a single pulsation frequency. For ground-based photometry, the 
aliasing problem and the higher noise level might also have con- 
tributed to the difficulty of detecting fx. 

We have found no evidence for modulation of the other fre- 
quencies, but we can comment on two other significant frequen- 
cies with amplitudes ~ 40ppm close to fx which are not in- 
strumental. These frequencies appear to be Fourier artifacts pro- 
duced by irregularities in the beating of fx and fx. As mentioned 
in SectionfJ) a multi-sine fit fails to converge if they are included 
in the solution. We have to mention here again that the broad 
filter band pass used by Most may be disadvantageous for such 
investigations. 

6.3. On the validity of interpolated frequencies 

To increase the resolution of our grid, we used linear inter- 
polation in temperature, luminosity, and polar magnetic field 
strength, respectively. It may be questioned if such an interpo- 
lation is a sensible approach. For a test we compared interpo- 
lated values with genuine models. For all genuine models m g 
enclosed by 2 additional genuine models m gt i and m Si x , we pro- 
duced frequencies at the position of m a by interpolation between 



m g i and m gi x- The deviation of the interpolated frequencies from 
their genuine values, normalized to the gap in log between 
the enclosing and the enclosed models, was calculated and over 
16000 frequencies were compared. Figure [9] shows a histogram 
of the results. It illustrates that our interpolation delivers a rea- 
sonable approximation to genuine models, as long as the gaps in 
between the reference models remain reasonably small. 

In the case of Grid 3 (see TableO we also repeated the fre- 
quency fitting procedure. We interpolated between 3a and 3c 
and omitted the 3b-models of Grid 3 . FigurefTOl shows the x 2 - 
statistics for our observed frequencies based on an interpolation 
in this much coarser grid (best fit with x 1 = 2.642). A com- 
parison of FigurefTUl with the upper panel of Fig. [6] again shows 
that interpolation is applicable for small gaps in between gen- 
uine models, which is also reflected in the results of the fit- 
ting procedure. Since not all modes can be approximated suf- 
ficiently through linear interpolation, as indicated by the outliers 
in Figure [9] it is necessary to check whether the modeled fre- 
quencies that fit the observations satisfy the linear assumption. 
Figure[JTJpresents the absolute (rather than the normalized) de- 
viation of the interpolated frequencies from the genuine frequen- 
cies for the 3/7-models. It shows that the interpolated values of 
all 5 modes that fit the observations are hardly deviating from 
their genuine values - the deviations lie well within the model 
uncertainties. 

It is important to note that all of the modes which fit the 
observed frequencies are actually not expected to be excited ac- 
cording to current model physics. This hints at remaining prob- 
lems concerning the stability analysis or the implemented driv- 
ing mechanism. 



6.4. Latitudinal amplitude dependence 

As discussed in ^5]the observed frequencies of y Equ are identi- 
fied as modes of 6 = 0, 1, 2, and 4. Those modes, however, have 
amplitude distribution on the stellar surface considerably devi- 
ating from that of a single Legendre function P( (cos 9) due to 
the strong magnetic effect (FigfT2l. The amplitude for/] (t = 1) 
and fn (( = 2) is more concentrated toward the polar regions than 
Pi (cos 9) and Px{cos9). It is interesting that the amplitude dis- 
tribution of fix {( — 2) on a hemisphere is not very different from 
that of f\ {( — 1) having very small amplitude near the equator, 
although f is antisymmetric and / symmetric to the equator. 

For fx, f 3 ({ = 4) and / ({ = 0) the amplitude distribu- 
tions are strongly concentrated around cos 9 ~ 0.65, drastically 
different from those for non-magnetic stars. Although the light 
variation from an I = 4 mode of a non-magnetic star is expected 
to suffer from a strong cancellation on the stellar disk, Fig. [12] 
indicates that our cases / and / seem hardly affected by this 
effect. 



7. Conclusion 

We have demonstrated that our models of roAp stars reproduce 
the observations in great detail. The temperature and luminos- 
ity parameters of our best fit are in goo d agreement with those 
expec ted from previous observations dKochukhov & Bagnulol 
12006b. although w e find a slightly higher mass. According to 
iRvabchikova et all (1 1997b the mean magnetic field modulus of 
y Equ is 4kG, which is half of the polar magnetic field 
strength of 8 kG as suggested by our best model. At present we 
can only speculate about a difference of internal magnetic field 
strengths, to which pulsation is sensitive, relative to surface mag- 
netic fields, which are accessible to spectroscopy. Furthermore, 
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Fig. 9. Upper panel: A histogram of the normalized deviation 
of 16050 frequencies, interpolated as a function of log T e g, from 
their values in genuine models. The bars show the number of fre- 
quencies deviating within a certain range. The dotted line gives 
the estimated upper limit of the model uncertainties. The dashed 
line shows the total percentage of frequencies, as indicated by 
the right-hand side ordinate, as a function of frequency devia- 
tion. 

Lower panel: A histogram of the absolute frequency deviation 
of 8523 frequencies, interpolated as a function of Bp, from their 
genuine counterparts. No normalization was necessary, since the 
interpolation was performed with a step size of 0.2 kG for each 
frequency. The dashed and dotted lines have the same meaning 
as in the upper panel. 
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Fig. 10. The result of fitting our observed frequencies to a model 
grid calculated by interpolation between 3a and 3c. Although the 
ife-models are omitted and no fit with^ 2 < 1 is found, the figure 
compares well to the upper panel of Fig. [6] Linear interpolation 
in Bp is therefore a reasonable approximation to genuine models 
for our fitting procedure. 
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Fig. 11. The absolute deviations of all interpolated frequencies 
at the position of 3b, using Grid models 3a and 3c, from the gen- 
uine frequency values. All 5 modes, which match the observed 
frequencies of y Equ best and are indicated by boxes, are devi- 
ating less than the model uncertainties. 



our models assume a simple, axisymmetric magnetic dipole, 
which may not sufficiently reflect a more complicated reality. 
The factor of 2 between (surface) magnetic field modulus and 
best fitting magnetic pulsation m odel appears again for the roAp 
star lOAql dHuber et al.l 120081) . another Most primary target. 
Certainly, more investigations of this sort have to be conducted 
in order to solve this problem. 

Pulsation amplitude changes forroAp stars were discussed in 
the literature as a consequence of limited mode life time or beat- 
ing frequencies. For y Equ we can clearly identify the beating 
frequencies, which seriously questions excitation mechanisms 
for roAp stars based on stochastic processes. 

Finally, we could identify the modes of all 7 frequencies de- 
tected so far in y Equ as I = 0, 1,2 and 4, with fa and fi being 
the harmonics of f\ produced by non-linear pulsation. Since we 



are not able to calculate non-axisymmetric modes, all of these 
frequencies are assumed to be (m = 0) p-modes. While we 
cannot rule out the possibility of excitation of frequencies with 
(m + 0), axisymmetric modes are most probable, because of the 
extremely long rotation period of y Equ and the assumption of a 
magnetic dipole. The assignment of any f-value to these modes, 
however, has to be understood as a convenient simplification. 
In the presence of strong magnetic fields a mode's oscillation 
behavior cannot be described by a single spherical harmonic. 
Consequently, we alert the reader that analysis techniques using 
this assumption, e.g. mode identification based on LPVs, should 
be adjusted for magnetic effects on the pulsation geometry. 
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Fig. 12. Latitudinal variation of the perturbation of radiative flux 
on the surface is shown for each mode matched to y Equ. (The 
case of fj is not shown because it is very close to that of /2O 
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poles). The plots are the real parts of the eigenfunctions for the 
flux perturbation at B P = 8 kG for model 3b (Table |3). Each 
curve is normalized as that the maximum is unity. 
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